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Nomenclature 


Subscripts 


c p specific heat at constant pressure, J/(kg-K) 

f v vapor friction factor 

h f 9 latent heat of vaporization, J/kg 

h 0 outer surface convective heat transfer 

coefficient, W/(m 2 K) 
k thermal conductivity, W/(m-K) 

Ls length of the liquid film to the point 

of complete condensation (Fig. 1), m 
I\ curvature, m' 1 

Ma w c / y/'yoRy^, Mach number 

jtif.m liquid mass fiow rate at z = 0, kg/s 
m-t total mass flow rate at z = 0, kg/s 

m mass flow rate, kg/s 

p pressure, Pa 

Q Q = 2 ^ r Jq q w dz— _ 

,TTl{)T( Cp t ( t in , 7Tl( t in ? T /',in ) i 3-Xial 

heat flow due to phase change, W 
q heat flux, W/m 2 

r radial coordinate, m 

R inner radius of the channel, m 

R y gas constant, J/(kg-K) 

Re 2w v (R - 6)p r /p v , axial vapor Reynolds 

number 

Rc ; 2v V '$ (R - 5)p v /p v , radial vapor Reynolds 

number 

T temperature, K 

v radial velocity, m/s 

it axial velocity, m/s 

Hr area-averaged axial velocity m/s 

z axial coordinate, m 


Greek Symbols 

n accommodation coefficient 

j3 momentum flux coefficient 

-yo c p /c t! , ratio of specific heats 

A dS /dz 

5 liquid film thickness, m 

p dynamic viscosity, Pa-s 

p density, kg/m 3 

<r surface tension, N/m 

ip inclination angle 


ent entrance 

in at z = 0 

i liquid 

men meniscus 

sat saturation 

v vapor 

w wall 

5 liquid film free surface 

crit at critical complete condensation 


Introduction 

Electronic cooling and other high heat flux applications re- 
quire a fundamental understanding of the condensation pro- 
cess in small diameter channels in order to optimize design 
configurations of devices used in these applications. Unlike 
conventional size passages in which surface tension effects are 
limited, surface tension in miniature size channels can have 
a significant role on the overall hydrodynamics and in par- 
ticular on the thin films that are believed to be the dominant 
mechanisms controlling the heat transfer characteristics. 

For forced convection condensation in miniature circular 
tubes the two phase flow regime maps for conventional sized 
tubes may not be relevant because of the role of surface tension 
in the hydrodynamics. Srinivasan and Shah (1997) point out 
that very limited data is available concerning the basic flow 
patterns for two phase flow with or without heal transfer in 
miniature circular tubes. 

Condensation in conventional size circular tubes comprise 
a well-defined body of work. Collier and Thome (1994) 
present a generally accepted description of the two phase 
flow patterns during forced convection condensation in con- 
ventional size horizontal tubes with co-current flow. The range 
of flow patterns present depend on the total energy convecled 
into the tube. In order of increasing velocity, the flow is found 
to be stratified, slug, plug or wavy, and annular. 

The objective of this work is to model annular film con- 
densation in miniature circular cylinder tubes where capillary 
phenomena can conceivably result in blocking of the tube cross 
section with liquid at some distance from the condenser en- 
trance. A physical description is illustrated in Fig. 1 (a). This 
phenomenon is referred to as complete condensation. An ex- 
periment with visualization was conducted and a mathematical 
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Figure 1: Description of the physical model for annular film 
condensation in a miniature tube, (a) Structure of two phase 
flow for complete condensation. (b)Coordinate system and 
conventions for film condensation model. 


model was developed on this subject in this work. 


Flow Visualization in Miniature 
Tubes 

A circular Pyrcx tube of 3.25 mm inside diameter and 5.00 mm 
outside diameter was used in a loop thermosyphon for flow 
visualization. The overall length of the horizontal test section 
was 302 mm with an insulated inlet section of 97 mm and a 
condenser length of 205 mm. Experiments were conducted 
over a range of heat loads from 4.9 to 1 0 W. 

Figure 2 is a photograph of the complete condensation phe- 
nomenon for the operating conditions of Q in = 7.7FT, T v = 
59.0°C. A liquid film on the lower wall is obvious closer to 
the location of complete condensation. A dramatic thinning of 
the film on the lower wall is also apparent. Immediately after 
the decrease in film thickness at the bottom wall, the liquid is 
seen to completely span the tube cross section. The closing 
off of the channel to vapor flow appears similar to a capillary 
tube meniscus. Note the slight inclination of the meniscus due 
to the effect of gravity. 


Mathematical Formulation 


pled vapor and liquid flows with shear stresses at the liquid 
free surface due to the vapor-liquid frictional interaction and 
surface tension gradient. The model is based on the following 
simplifying assumptions: 

1 . The vapor is saturated and there is no temperature gra- 
dient in vapor in radial direction. 

2. Heat transport in the thin liquid films is only due to 
conduction in the radial direction. 

3. Inertia terms can be neglected for the viscous flow in 
the liquid films with low Reynolds numbers. 

4. Force on liquid due to surface tension is much greater 
than the gravitational force and therefore the liquid is 
distributed onto the walls in a film of locally uniform 
thickness. 

5. The solid tube wall is infinitely thin so that its thermal 
resistance in the radial direction can be neglected as 
well as the axial heat conduction. 

The cylindrical coordinate system used is shown in Fig. 
1(b). Both the vapor and liquid flow along the c -coordinate. 
The physical situation should be described taking into consid- 
eration the vapor compressibility and the vapor temperature 
variation along the channel. Also the second principal radius 
of curvature of the liquid-vapor interface should be accounted 
for (in the equation relating vapor and liquid pressures) while 
it is usually neglected in modeling of film flows in tubes of 
larger diameters. The mass and energy balances for the liquid 
film shown in Fig. 2(a) yield: 



T w - T s 
In [R/(R-S)\ 



7h ( T() 


( 1 ) 


where T ( is the area-averaged liquid temperature for a given 
z location. For consideration of subcooling in the condensed 
liquid, T ( is found from an average area given by 


1 [R 2 -{R- <S) J ] 


(2) 


where T t {r) is the assumed liquid film temperature profile 
given by the temperature distribution in a cylindrical wall. 


A summary of the analysis is given below. For the complete 
development of the mathematical model see Begg et al . (1998). 

A steady state mathematical model of condensation which 
leads to complete condensation is presented and includes cou- 


T ( (r) = T 6 + 


Tw - Ts 


In 


R — 5 


The derivative of T t is approximated by 


( 3 ) 
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Figure 2: Photograph of annular film condensation resulting 
in complete condensation in a 3.25 mm inside diameter tube. 
Q itl = 7.7JF, T v = 59.0°C 


(1T( 

dz 


dT s ( dT w _ dTjj > 
+ \ dz dz ) 
'R-5' 


•-(nr) 


[fl 2 - (fl - <5) 2 ] 

2 («-*)*! 


(4) 


The momentum conservation for viscous flow in a liquid 
film in which the inertia terms are assumed to be negligible is 


if (r2p) 




I — + pegsuvp 


) 


(5) 


The boundary conditions for the last equation are the non-slip 
condition at r — i? and specified shear stresses at the liquid- 
vapor interface due to the frictional liquid-vapor interaction, 
r ltVt and the surface tension gradient related to the interfacial 
temperature gradient along the channel. 


uv| r =/? = 0 


( 6 ) 


dun . 


/*<■ 


da dTs' 
r T( ' v ~ dT~dz. 


= E 


(7) 


where T 6 is the local liquid-vapor interface temperature and 
^ is the Marangoni effect. 

Taking into account the effect of the condensation process 
on the shear stress term, an expression from the Munoz-Cobol 
el al. (1996) analysis is used for t ( , v . The equation for the 
axial pressure gradient in the liquid is 


dpc 

dz 




+pc 


1 

2np( 


^--Th t , in ) + E(R - 5)F 

9 


[ f? 4 {R — S) 2 f (R-S) 1 

16 + V 8 4 ) 


( 8 ) 


where 


F = 


(R-Sl 


(, R m — 
V i?-5 + 2/ 4 


(9) 


The pressure difference between the vapor and liquid 
phases is due to capillary effects (Faghri, 1995). 



Tzs m ( aaa h)}- pd (10) 

The term with cosine in the right-hand side of this equation is 
due to the second principal radius of the interfacial curvature 
for a cylindrical film. Introducing an additional variable 


dS/dz = A (11) 

eq. (10) can be rewritten as follows. 


dA 

dz 





Pv - p{ + Pd 
a 


cos(atanA) N \ 
R-S ) 


( 12 ) 


The compressible quasi-one-dimensional momentum 
equation for the vapor flow in the form suggested by Bankston 
and Smith (1973) is modified to account for non-uniformity 
of the vapor cross-sectional area of the liquid-vapor interface 
following Faghri (1995). 


dpr 

-j- = p v g snv+ 
dz 

-J- [4 (-A P.^-4,) - frPrw'MR ~ 6) + 

-H v L dz 

2n{R - 6)p v vlj sin (atanA)] (13) 

with p v = 1.33 for small radial Reynolds numbers. 

The perfect gas law is employed to account for the com- 
pressibility of the vapor 


Therefore 


Pv 

RyTv 


(14) 


dp v _ 1 (dp v 1 p v dT v \ 
1h-W y \dz T v n dz) 1 j 

The saturated vapor temperature and pressure are related by 
the Clausius-Clapeyron equation which can be written in the 
following form. 
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dT v __ dp v RgT v ^i^) 

dz dz p v hjg 

The seven first order differential equations Eqs. (1), (8), 
(11), (12), (13), (15), and (16) include the following seven 
variables: 5, A, pf, Q, pv, Pvi and T v . Therefore, seven 
boundary conditions are set forth at z = 0 


6 = Sin 

( 17 ) 

> 

It 

O 

(18) 

2 a 

in R-S in +Pd 

(19) 

Q = 0 

(20) 

t .in “ Pv, sat(7r, in) 

(21) 

Pi', in 

in t-> rp 

t\y j- v,in 

(22) 

pt 

II 

(23) 


In boundary condition (17), 6, n is defined from the condition 
in the adiabatic zone at the entrance to the condenser where the 
liquid and vapor pressure gradients are equal. This condition 
is satisfied by solving equations (8) and (13) for Sin for the 
case of Q =0, v Vt s =0, dA v /dz =0 and dTs/dz =0. There are 
also parameters mc, t n and w v ^ n and an additional variable 
T 6 involved in this problem. They will be considered using 
additional algebraic equations. The parameter should 

be found using a constitutive condition at the entrance of the 
condenser. 


Ttiejn = rht - Qt/hfg (24) 

where Qt is the total heat load of the condenser. Also w v ,in — 
Qt/(hfgPvA Vt in)- 

The liquid- vapor interface temperature, Ts, differs from 
the saturated bulk vapor temperature because of the interfacial 
resistance and effects of curvature on saturation pressure over 
liquid films. The interfacial resistance, is defined as (Faghri, 
1995): 


The following two algebraic equations should be solved to de- 
termine Ts for every point along the 2 -direction. The relation 
between the saturation vapor pressure over the thin evaporating 
film, (p S at)<5, affected by the surface tension, and the normal 
saturation pressure corresponding to Ts, p^i(Ts) f is given by 
the extended Kelvin equation (Faghri, 1995): 


(Psai)t5 — Psai{Td) 

(psat )<S ~ psaijTs) ~ &R 


exp 


-Pd 


p(RgTs 


(26) 


where K is the local curvature of the liquid-vapor interface 
defined by the term in outer brackets in eq. (10). Noticing 
that under steady stale conditions, qs is due to heat conduction 
through the liquid film, it follows from this condition and eq. 
(25): 


Ts — T w + 


(R-S) 

k( 


In 


R 

R — 6 


hfy Pv 
s/lnRg L V% 



(27) 


Equations (26) and (27) determine the interfacial temperature, 
Ts , and pressure, (p sa t)«$, for a given vapor pressure, p v = 
p r sal (T t .), temperature of the solid-liquid interface, T u ., and 
the liquid film thickness, S. 

For the case of variable wall temperature, an additional 
variable, T w , and first order ordinary differential equation must 
be added to the seven previously specified. 

T w is local wall temperature and can vary along the con- 
denser length. If the convective heat transfer coefficient at the 
outer tube wall, h a and the cooling liquid temperature, T 0 c , 
arc known, the local wall temperature can be defined using an 
energy balance. 


h 0 {T U ’ — Too) — 


1 dQ 
2 ttR dz 


(28) 


Another boundary condition is also required and that is 
given by 


2 =0, T u . = T u ., in (29) 

Numerical Procedure 


/ 2a \ hjg 

Pv 

(Psal)rf 

\2 — a) y/lnRy 


JTs . 


where p v and {pm)* are the saturation pressures correspond- 
ing to T r and at the thin liquid film interface, respectively. 


Equations (1), (8), (13), (15), and (16) with corresponding 
boundary conditions have been solved using the standard 
Runge-Kutta procedure. Algebraic equations (26) and (27) 
with two unknowns, (psn^ and Ts , have been solved numeri- 
cally for every point on z using a standard numerical procedure 
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(Wegstein’s iteration method). All the unknown variables were 
found with the accuracy of 0.0005%. During the numerical 
procedure the interval 0 < z < L$ was divided into at least 
500 parts, and the thermophysical properties of the saturated 
vapor and liquid were recalculated for each of the parts at the 
corresponding vapor temperature T v (z). 


Results and Discussion 


Figure 3 (a) shows variation of the liquid film thickness along 
the condenser with a constant wall temperature of 340 K. The 
vapor (TV, = 363 K) pressure drop over the condensation 
length is insignificant compared to the liquid pressure drop as 
shown in Figs. 3 (b, c). Condensation is more intensive in the 
region where the film thickness is at the minimum. For the 
cases of Q in = 8 W and Q in = 10 W, representing complete 
condensation and critical complete condensation, respectively, 
the incoming vapor is totally condensed. This is not true for 
Q in - 12 W, the incomplete condensation case. Figure 3 (d) 
shows the cumulative sum of the heat removed from the vapor 
by condensation. For Q in = 12 W, not all ol the heat load 
coming into the condenser is rejected at the termination of the 
calculation. 

Results are obtained for the convective cooling boundary 
condition shown in Fig. 4. The wall temperature is non- 
uniform with a significant increase just prior to the point of 
complete condensation, coincident with the thinning of the 
liquid film, as shown in Fig. 4 (d). Decrease of the convective 
heat transfer coefficient, h 0 > from 20,000 to 15,000 W/m 2 K 
resulted in an increase of the wall temperature that shortened 
the condensation length, as shown in Fig. 4 (a, d). 

Summary of the Complete Condensation Phenomenon 

Features of complete condensation in the annular film regime 
include the segregation of liquid and vapor into di stinct regions 
within the tube. Initially, the annular film thickness increases 
gradually in the downstream direction. After reaching the 
local maximum value, it suddenly decreases to its minimum 
thickness and then dramatically increases to the radius of the 
tube. The region at the convergence of the liquid film is similar 
to the classic capillary meniscus in appearance. This forms a 
well defined location marking the transition to single-phase 
liquid flow. Downstream of the location of the meniscus- 
like surface, the entire cross section of the tube is occupied 
by liquid. Any additional heat removal downstream of this 
location only results in subcooling of the liquid. The axial 
distance from the inlet of the condenser tube at z. = 0 to the 
location at which the entire cross section is filled with liquid 
is defined as L*, the condensation length. 



z(jim) 


Figure 3: Annular film condensation in circular tube with con- 
stant wall temperature. R = 1.5mm, mt = O.Olg/s, TV- = 
340 K. TV = 363 A' (a) film thickness (b) vapor pressure (c) 
liquid pressure (d) cumulative heat rate rejected 
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Conclusions 


Based on the numerical results and visual observations, the 
following conclusions have been made. 


1 . Length of film condensation in miniature tubes is very 
restricted due to surface tension effects that result in 
complete condensation a short distance from the con- 
denser inlet. 




z (pm) 


Figure 4: Annular film condensation in circular tube with 
convective boundary condition. R = 1.5mm, m* — 
0.01tf/.s, T v%in = 323 Jv, 7k = 300/V, Q in = 5.83TT 
(a) film thickness (b) vapor pressure (c) liquid pressure (d) 
wall temperature 


2. The complete condensation length is a non-linear func- 
tion of the vapor- wall temperature difference, increas- 
ing with the temperature difference approaching zero. 

3. Pressure drop in the vapor phase was usually insignif- 
icant compared to the pressure drop in liquid. Con- 
sequently, the vapor temperature variation along the 
condenser was infinitesimal. 

4. Variation of the wall temperature can be significant for 
convective cooling boundary condition. 

5. Complete condensation in small diameter tubes is a 
truly steady phenomenon. 

A very distinct transition is anticipated between the steady 
phenomenon of complete condensation and the next two-phase 
flow regime, which is referred to here as incomplete conden- 
sation. Tt is anticipated that incomplete condensation will be 
characterized by fluctuations of the liquid- vapor interface in the 
annular film which result in bubbles of vapor passing down- 
stream of confluence of the liquid film. The meniscus-likc 
feature will no longer be stationary and an unsteady two-phase 
flow bubbly-type flow will exist in the tube. 
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